Coherence of a qubit stored in Zeeman levels of a single optically trapped atom 



O 



5^5 



Wenjamin Rosenfeld,''^'Lj Jiirgen Volz,'H| Markus Weber/ and Harald Weinfurter''^ 

Fakultdt fiir Physik, Ludwig-Maximilians-Universitdt Munchen, D-80799 MUnchen, Germany 
Max-Planck-Institut fur Quantenoptik, D-85748 Garching, Germany 

We experimentally investigate the coherence properties of a qubit stored in the Zeeman substates of the 
1 M, F = 1 hyperfine ground level of a single optically trapped Rb atom. Larmor precession of a sin- 
gle atomic spin- 1 system is observed by preparing the atom in a defined initial spin-state and then measuring 
the resulting state after a programmable period of free evolution. Additionally, by performing quantum state 
tomography, maximum knowledge about the spin coherence is gathered. By using an active magnetic field sta- 
bilization and without application of a magnetic guiding field we achieve transverse and longitudinal dephasing 
times of T^ = 75..150/XS and Ti > 0.5ms respectively. We derive the light-shift distribution of a single atom 
in the approximately harmonic potential of a dipole trap and show that the measured atomic spin coherence is 
limited mainly by residual position- and state-dependent effects in the optical trapping potential. The improved 
understanding enables longer coherence times, an important prerequisite for future applications in long-distance 
quantum communication and computation with atoms in optical lattices or for a loophole-free test of Bell's 
inequality. 
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I. INTRODUCTION 

Quantum memories for the storage and retrieval of quantum 
information play an outstanding role in future applications of 
quantum communication, such as quantum networks and the 
quantum repeater |[l|] . There, ground states of trapped atoms 
or ions, are ideal candidates, as the interaction with the en- 
vironment is weak and can be controlled with high accuracy. 
Although in such systems coherence times of the order of sev- 
eral seconds have been observed [2^-^, storage and retrieval 
of single quantum excitations was shown to reach maximum 
times of only few 100/xs to ms i5|-|9|], in which storage of 
a complete polarization qubit state represents a greater chal- 
lenge. In order to further prolong the quantum storage time 
a detailed understanding of dephasing and decoherence pro- 
cesses is indispensable. 

Optical dipole traps lUOll nowadays are a well-established 
tool for the controlled manipulation of internal and external 
quantum states of neutral atoms 1 3l Ift Isl [lll - [l6ll . Such traps 
provide almost ideal conservative trapping potentials com- 
bined with low heating rates, resulting in long atomic coher- 
ence times. However, two kind of effects significantly limit 
the achievable coherence time: (i) Off-resonant spontaneous 
Raman scattering from the dipole laser beam entangles the 
qubit of freedom with some degree of freedom of a single 
scattered photon II17I - U9I1 . This kind of light-matter interac- 
tion leads to decoherence in the most general sense as the sys- 
tem under investigation (atomic memory qubit) gets entangled 
with the environment (scattered photon), (ii) In addition, if the 
trapped atom is not in the vibrational ground state, the ther- 
mal motion together with residual state-dependent effects of 
the optical trapping potential will lead to dephasing of an ini- 
tial atomic spin state. Although the temporal spin-evolution 
in an individual experimental realization is strictly coherent. 



the observed ensemble average over many experimental runs 
may have a large scatter due to different initial conditions. 



In our experiment quantum information is stored in the Zee- 
man sublevels of the S^Sin, F = 1 hyperfine ground level of 
a single '^^Rb atom, localized in an optical dipole trap ll20ll . 
For a variety of applications in long-distance quantum com- 
munication, e.g., the generation of long-distance atom-photon 
Oa llMl and atom-atom entanglement, as well as the closely 
related task of the remote preparation of an atomic quantum 
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memory 112 111 , only a spin-1/2 subspace of the 5 Sin, F ~ 1 
hyperfine ground level is addressed. More precisely, the 
atomic memory qubit is encoded in the \F = I, mp = — 1) and 
\F — I, nif = +1) Zeeman sublevels. The remaining third 
sublevel \F =\, nif = 0) is not directly used for qubit storage, 
however, the coherent evolution of the total angular momen- 
tum F = 1 in a magnetic field (Larmor-precession) can lead 
to its population, thereby reducing the fidelity of the stored 
state. In order to extract information how a stored quantum 
state becomes mixed and also to distinguish coherent from in- 
coherent processes, the temporal evolution of the full spin-1 
density matrix of the S'^Sin, F = I hyperfine ground level has 
to be investigated. 



In this contribution we analyze decoherence and dephasing 
mechanisms and their relevance for the storage of quantum 
information in single optically trapped ^^Rb atoms. In detail, 
in Sec. |II]we develop a model accounting for state-dependent 
effects of the optical trapping potential. In order to achieve 
spin-coherence times of several 100 /xs we implement an ac- 
tive magnetic field stabilization (see Sec. HIH) . Applying par- 
tial quantum state tomography of the total angular momentum 
state 5^ Sin, F —Iwe then investigate in detail dephasing and 
decoherence mechanisms. Finally, in Sec. |IV]we summarize 
our major findings and give an outlook how longer spin co- 
herence times could be reached. 



II. THEORETICAL ANALYSIS OF DECOHERENCE AND 
DEPHASING MECHANISMS 

A. Decoherence due to spontaneous Raman scattering 

The most influential scattering process occurs due to inter- 
action of the atom with the light of the dipole trap. In our case 
the dipole trap is generated by a single, sharply focused laser 
beam at a wavelength of A = 856 nm. Despite of the large de- 
tuning to the first dipole allowed transitions S^Sin -> S^Pj /2 



(A = 795 nm) and 5^Si/2 ^ 5^^3/2 (A = 780nm) in **^Rb, 
there is still a finite probability to spontaneously scatter light 
from the dipole laser beam. This scattering process consists 
of two important parts. Elastic (Rayleigh) scattering occurs 
when the atom returns to the same state after emission of a 
photon. In our case this happens at a rate of 17.7 Hz ll22ll . 
When the final and initial states are different, even though the 
states might be degenerate in energy, this process is called 
spontaneous Raman scattering. As demonstrated previously 
1 17-19], it is this scattering process which entangles, e.g., the 
polarization and/or the frequency of a single scattered photon 
with the internal spin state of one or many atoms. Obviously, 
only spontaneous Raman scattering will lead to spin relax 
ation. Far from the atomic resonances 5^S 



>l/2 



5^Pi/2 and 



5 Si/2 -^ 5 P3/2 spontaneous Raman scattering is strongly 
suppressed due to destructive interference of amplitudes in the 
different excitation and decay channels ||23[1 . Its rate is given 
by 
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The optical dipole trap is based on a spatially varying, light- 
induced energy shift of the atomic ground levels. Assuming 
that all relevant detunings are larger than the hyperfine ground 
state splitting we obtain the light-shift for the Zeeman sub- 
levels \F — 1, nif) of the S^Sin ground level of (see e.g. 1101 ) 
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For linear polarization of the trap light (j3^ ^Q) this energy 
shift is equal for all magnetic sublevels ntf of both hyperfine 
ground levels 5^S\|2^F = 1 and S^Sin^F = 2 and therefore 
ideally suited as a state-insensitive trapping potential. How- 
ever, in the case of circular polarization (.^ = +1 for (7+ and 
3^ = —\ for (7^) the shift becomes state-dependent, lifting the 
degeneracy of the Zeeman sublevels. This additional effect on 
the magnetic sublevels nif ^ for circularly polarized light is 
formally equivalent to a magnetic field Be pointing along the 
propagation direction z of the dipole laser beam, and is called 
vector light shift (also Zeeman light shift) |24]. In contrast to 
an external magnetic field which can be considered homoge- 
neous over the microscopic volume of the optical dipole trap, 
the vector light shift depends on the intensity and thus on the 
position of the atom in the trap. For a typical temperature 
of the atom of 100..150;UK the thermal motion will lead to 
a non-negligible variation of light shifts, which significantly 
influence the dephasing of stored quantum information. This 
dephasing mechanism will be analyzed in the following. 



where © is the angular frequency of the dipole laser, / its in- 
tensity, Aoi and A£)2 its detuning with respect to transitions 
to 5^Pi/2 and 5^^3/2 levels. Yd is the total resonant scatter 
rate of the respective D\ and D2 transition, and (Od the an- 
gular frequency of the corresponding atomic resonance. The 
rate for spontaneous Raman scattering via higher-lying rp-Pin 
and rp-P^n levels (« > 5) is negligible due to the large detun- 
ing of the dipole laser For our typical dipole trap parameters 
(A — 856 nm, 7=1 .56 • 10^ W/m^) the incoherent scatter rate 
is 0.11 Hz. Thus, the atomic spin will relax on a time scale 
of about 10 s, which is comparable to the lifetime of captured 
atoms in the trap 120(1 . We conclude that in for our exper- 
iments spin relaxation due to spontaneous Raman scattering 
does not limit the coherence time. 



1. Coherent state evolution 

As a first step we calculate the coherent temporal evolution 
of the spin-1 state |^(f )) in a constant effective magnetic field 
Beff- The interaction Hamiltonian Hefj is given by 
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where F is the operator of the corresponding angular momen- 
tum F — I, Hb = 'h-2n ■ 1.4MHz/G is the Bohr magneton 
and gf = —1/2 the Lande factor The effective magnetic field 
Bgff is given by 



B. Dephasing of the atomic spin 

In contrast to the possibility that the atomic spin decoheres 
due to entanglement with the environment, the coherent in- 
teraction with fluctuating external magnetic fields B leads to 
dephasing of stored quantum information. Additional spin- 
dephasing for trapped atoms can also be caused by the ther- 
mal motion of the atom in a state-dependent trapping potential 
(see Sec. HTbH i. 
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with /(J the intensity of the circularly polarized component 
of the dipole trap beam. For convenience we set Beff — 
Beff{bxSx + byBy + b,e^), where b^ = ^/l^l}^cos{(j)), by — 
VT^sin(0). 
For this field configuration we obtain the eigenstates 
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with corresponding eigenvalues +hcOL, 0, —hcoi, and Lar- 
mor frequency col = jrjXBgFBeff- These are mp = +1,0,-1 
eigenstates with respect to the direction of the effective mag- 
netic field. Finally, for the time evolution of an arbitrary 
5^Si/2, F ~ 1 state we obtain 
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where the amplitudes are given by c±\ = (<I>±i|^(f = 0)), 
Co ~ ('I>o|1'(f = 0)). Obviously, the Larmor precession in an 
effective magnetic field Bf.ff is thus not necessarily limited to 
the qubit space {|1,-1), | lj+l>}. 



Il22ll . This can be understood as the atomic spin, which rotates 
at a finite Larmor frequency, can not follow the increasing fre- 
quency of the field fluctuations which therefore average out. 
Thus, as the magnitude of rapid fluctuations in our experiment 
is small (< 0.3 mG) and due to the additional l/Q? suppres- 
sion, we conclude that the influence of rapid oscillations is 
negligible, particularly when compared to the effect due to 
the slowly varying field component. 

For fluctuations of the magnetic field which are slow on the 
time scale of a single experimental run, the field can be con- 
sidered constant and the atomic state will evolve according to 
Eq. (|6]l. However, the field can vary between repeated exper- 
imental runs. This inevitably leads to different evolutions of 
the atomic state and therefore the observed average state pop- 
ulations are washed out. This dephasing can be modeled by 
first calculating the temporal evolution of the considered state 
I'P) in a constant effective magnetic field Beff according to 
Eq. (|6]), and incoherently averaging over states resulting from 
the distribution of different magnetic fields corresponding to 
different experimental runs. 

In our spin-precession experiments we start with an initial 
state 1^) and let it precess for a time t giving |^) (f ). Then the 
population of a chosen analysis state |^a) is obtained from 
the overlap with the precessed state, averaged over all possible 
evolutions. It is given by 



2. Fluctuations of the magnetic fields 

As a first step in the analysis of dephasing mechanisms one 
has to know the magnitude of magnetic field fluctuations at 
the relevant time scales. Typically, our experiments require 
preservation of the atomic quantum state for several microsec- 
onds 1 25]. In this work we consider times up to 200 /is, defin- 
ing two important frequency ranges. The first range contains 
frequencies, where the magnetic field varies rapidly on the ex- 
perimental time scale (i2 > 2;r • 2.5kHz). In the second range 
we have £1 <2n ■ 2.5kHz, i.e., the field can be considered 
constant within one experimental run, but will vary between 
different runs. With the help of a magneto-resistive sensor 
(accessible frequency range: DC — 60 kHz) we were able to 
quantify the magnitude of magnetic fields at different frequen- 
cies, see Sec. IIIIBI We found that the strongest fluctuations 
were at low frequencies (< 200 Hz) while faster fluctuations 
were relatively small (< 0.3 mG rms within the bandwidth of 
60kHz). Magnetic field fluctuations at different frequencies 
affect the atomic state in different ways as will be discussed 
below. 

In the case where the field fluctuates rapidly on the exper- 
imental time scale (and also rapidly compared to the Larmor 
frequency O)/,), the evolution of the atomic state will follow the 
average field B with small oscillations around the main trajec- 
tory. If the average field B is constant, only those deviations 
will lead to dephasing. The magnitude of these deviations 
can be estimated by solving the time-dependent Schrodinger 
equation in a field modulated at a frequency £1^ (Oi and was 
found to drop with increasing modulation frequency as 1 /Q? 
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where the Pj{Bj) are the normalized distributions for the /' = 
x,y,z components of the effective magnetic field. If, e.g., the 
average values are B^ = By = B^ = 0, and the z-component of 
the effective magnetic field follows a Gaussian distribution 
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then, after initially preparing the atomic spin state Iv'^i) = 
-j=(|l,-l) +e''''|l,+l)), we get for the probability to stay 
in this state after a time t 
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Here T^* 



can be associated with the transverse co- 
. Note that the decay 
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herence time of two-level systems (2 
for this noise model dHJ is not exponential 

As a second example we consider the dephasing of the spin- 
states 1 1/2) = |1,±1) in a fluctuating magnetic field along the 
X-axis. By averaging over the Gaussian distribution px{Bx) as 
in ([8]l we find 
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Here, the state population decays with a time constant T = 
, which is twice as long as 7^2*. P{\f/2){t) approaches 
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Figure 1: Probability density p{ABa) of the optically induced mag- 
netic field in a three-dimensional harmonic trap resulting from ther- 
mal motion of the atom. The distribution is plotted as a function of 
the deviation AB^ = B^ — B^ of the field B^ from its maximal value 
B^y at the bottom of the trap. The curve was calculated according to 
Eq. l f76t . assuming a trap depth Uq = kg ■ 650 /iK, average atomic 
temperature T = 150/xK and 1% fraction of circularly polarized trap 
light. 



the limit of |, as the corresponding spin evolution leaves the 
qubit subspace {|1, — 1),|1,+1)}. Such a situation does not 
occur in two-level systems and represents a more complex de- 
phasing scenario. For more general cases the integral in (|7]i is 
not analytic. 



3. Effective fluctuation of the light shift 

The second part of the effective magnetic field - the vector 
light shift - results from the circularly polarized component of 
the dipole trap light and is proportional to its intensity. Due to 
the thermal motion of the atom in the trap, in each realization 
of the experiment the atom will be found at a random posi- 
tion within the trapping potential and will therefore be subject 
to a different light shift (see inset of Fig. [1]). Here we shall 
consider the case, where the atom can be considered static 
within one experimental run. In our experiment the shortest 
oscillation period in the trap is 44 /is, therefore this assump- 
tion is strictly valid only for shorter time scales. Nevertheless, 
even for longer experimental times it represents a worst-case 
assumption, since a static field which changes from experi- 
ment to experiment leads to a stronger dephasing than a field 
of the same amplitude which fluctuates on the experimental 
time scale or faster (see Sec. IIIB2I I. 

The distribution of positions depends on the thermal energy 
of the trapped atom and maps directly onto a distribution of 
the induced magnetic field B<j. To derive this distribution one 
has to know the 3D distribution of the corresponding poten- 
tial energy. For a thermal energy sufficiently lower than the 
trap depth the potential can be considered harmonic and the 
distribution p{AB(j) is calculated as follows. 

The potential energy U of a 1 -dimensional harmonic oscil- 
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where E is the fixed total energy of the motion and (p is the 
phase of the oscillation. We define the potential energy U be- 
ing non-negative, with U — Q at the bottom of the harmonic 
trap. If at a certain random moment of time the potential en- 
ergy is measured, some random realization for the phase (p £ 
[0..2;r] will be found, where every value of (p has equal prob- 
ability. Then the probability to obtain a value of the potential 
energy within the interval [U, U + dU] is pE{U)dU <^ d(p{U). 
Here, pe{U) is the probability density for a given total energy 
E, given by 
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In thermal equilibrium the total energy E follows a Boltzmann 
distribution. The corresponding thermal distribution of the po- 
tential energy U is given by 

pw{U)= [ dEpE{U)-^expi--^) 
Jo ksT kBi 

For three dimensions the thermal distribution of U is obtained 
from a convolution of the three independent ID distributions 
II22I1 . resulting in 
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The corresponding kinetic energy Eiii„ follows the same dis- 
tribution. 

Here, it is worth mentioning that the distribution ( fT4] i differs 
from the well-known Maxwell-Boltzmann distribution, which 
is pmb{E) = j^j^E^expi-jMj,). According to the Virial 

theorem the average of the potential energy is half of the total 
energy: {U) — jE. This relation might suggest that the poten- 
tial and kinetic energy follow the same distributions as the to- 
tal energy E. However, that is not the case. The Virial theorem 
makes only a statement about average values, while the con- 
sidered distribution describes the probability to find a certain 
value of potential energy at a randomly chosen point in time 
(and is therefore not ergodic). Instead, it can be easily verified 
that the convolution of the distribution (fl4l i of the potential 
energy U and of the identical distribution PioiEkin) of the ki- 
netic energy Ei^j„ gives the expected Maxwell-Boltzmann dis- 
tribution of the total energy 

/ PioiU =E- Eki„)pj,D{Ekin)dEki„ = pmb{E). (15) 

Now we are able to derive the distribution p{AB(j,T) of 
the optically induced magnetic field Ba for a thermal atom 
at a given temperature T. For convenience we introduce the 
maximal value of the optically induced magnetic field at the 



bottom of the trap Z?^ according to Eq. (|4|i. For our typi- 
cal trap parameters 1 % of circular admixture in the polariza- 
tion of the trapping light results in B% ~ lOmG. The rela- 
tion between the induced magnetic field and the potential en- 
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Finally, we define the deviation 



AB(j = B^—B(j from the maximal value B^ at the trap center 
(Fig. [T] inset). Using these relations we obtain the distribu- 
tion of the optically induced magnetic field for an atom at a 
temperature T 
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This distribution is shown in Fig. [T]for typical experimental 
parameters. 

Expression (fTSI l directly relates the thermal motion of the 
trapped atom and the circular admixture in the polarization 
of the trapping light to a fluctuating effective magnetic field. 
These fluctuations can have a serious impact on the achievable 
coherence times as will be shown in Sec. |III] 



III. EXPERIMENTAL ANALYSIS OF DEPHASING 
MECHANISMS 

A. Single atom trap 

In our ex peri ment a single ^^Rb atom is stored in the optical 
dipole trap Ill9ll20ll . which is loaded from a laser-cooled cloud 
of 10^.. 10^ atoms of a shallow magneto-optical trap (MOT). 
The conservative optical trapping potential is created by a 
focused Gaussian laser beam (waist wq = 3.5 /xm; Rayleigh 
range zr = 4-5 fim) at a wavelength of 856 nm, thereby de- 
tuned far to the red of any atomic transition from the atomic 
ground level. For a typical power of P = 30 mW we achieve 
a potential depth of Uq ~ ks ■ 650/iK, corresponding to ra- 
dial and axial trap frequencies of (Or — 27r- 22.7kHz and 
©J = 2n ■ 1.25 kHz, respectively. This trap provides a stor- 
age time of several seconds, mainly limited by collisions with 
the thermal background gas ll20ll . The fluorescence light of the 
trapped atom is collected in a confocal arrangement by an ob- 
jective, coupled into a single mode optical fiber and guided to 
two single photon counting avalanche photo-detectors (APDs) 
allowing also a polarization analysis of single photons. The 
presence of a single trapped atom is inferred by detecting flu- 
orescence light. The bare detection efficiency for a single pho- 
ton emitted by the atom is 2 • 10^^, including coupling losses 
into the single-mode optical fiber and the limited quantum ef- 
ficiency of the single photon detectors. 



B. Active magnetic field control 

A crucial requirement for achieving long atomic coherence 
times is a precise control of the magnetic field in the region 
of the optical dipole trap. For this purpose the fields in our 
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Figure 2: (Color online) Schematic of the experimental setup (not 
to scale). The dipole trap beam is focused into a vacuum cell by 
means of a high NA objective. The same objective collects the light 
emitted by the atom into a single mode optical fiber. The magnetic 
field sensor positioned on the cell surface is used to give a feedback 
signal to compensation coils for field stabilization. For simplicity, 
only one pair of compensation coils is shown. 



experiment are actively stabilized. The magnetic field is con- 
tinuously monitored by a three-axis magneto-resistive sensor 
(Honeywell HMC1053), which is located outside the vacuum 
glass cell at a distance of 25 mm from the position of the 
trapped atom (see Fig. |2l). On a short time scale the pre- 
cision of this sensor is limited by electronic noise (typically 
< 0.1 mG rms within the effective bandwidth of 60kHz). A 
more significant problem is saturation of the sensor by the 
strong fields during loading of the MOT. In this case a re- 
magnetization of the sensor is required which limits the pre- 
cision to 0.5..1mG on long time scales. Using this sensor 
we have measured the fluctuations of the magnetic fields and 
could identify two main sources. The largest part of the fluc- 
tuations is due to currents drawn by the Munich underground 
train line passing at a distance of about 60 m from our lab- 
oratory. The time scale of these fluctuations is 30s..lmin 
with a peak-to-peak amplitude of 20..25mG on the vertical 
and 6..8mG on the horizontal axis. The second major contri- 
bution arises from the 50 Hz mains supply producing fluctu- 
ations of about 1..2mG peak to peak on each axis. For fre- 
quencies higher than the fourth harmonic of the power line 
frequency (200 Hz) the fluctuations were found to be on the 
order of < 0.3 mG rms. 



The signal from the magnetic field sensor is fed back 
to compensation coils by means of a servo loop. The in- 
tegration time constant was set such that an active band- 
width of about 200 Hz was reached, sufficient to suppress 
the effects of underground trains and the power supply line. 
Given the fluctuations of external fields described above, our 
active magnetic field stabilization achieves an rms stabil- 
ity of (0.92, 0.77, 0.83) mG for the three axes, including re- 
magnetization precision of the sensor, crosstalk between dif- 
ferent axes (< 3.5%) and magnetic field gradients between 
the position of the trapped atom and the position of the sensor 
(< 5%). 



C. State preparation and detection 

To study dephasing of a single atomic spin 5^51/2, F — 1 
we initialize the atom with high fidelity in a well-defined 
state of our choice. This is realized by first entangling the 
atomic Zeeman-states \F ~ \,mp = — 1) and \F = I, nip — 
+ 1) with the polarization states |ff+) and |a^) of a sin- 
gle emitted photon IQ, [iM |2lll - the entangled state reads 
!>!'+) = -i=(|l,-l)|a+) + |l, -hi) |(T"))- and then projecting 
the atom onto the desired spin-state via a polarization mea- 
surement of the photon. Thus, a measurement of the pho- 
ton in, e.g., the cr*-basis (circularly polarized) leaves the 
atom in the |l,Tl) state and a measurement in ///V-basis 
(horizontally/vertically polarized) projects the atom into the 
-^(|1,-1)±|1,+1)) states, respectively. 

After preparation, the atomic spin freely evolves for a de- 
fined period of time in the applied magnetic field B. During 
this time all lasers except for the one used for the dipole trap 
are switched off and the magnetic field is stabilized to a pre- 
selected value. Finally, after a given time, the atomic state 
detection procedure 111911 is applied, allowing us to determine 
the projection of the atomic state on any superposition in the 
{|1, — 1),|1,+1)} subspace. By these means we can directly 
observe the temporal evolution of selected atomic states (see 
Fig.0. 



D. Analysis of the state evolution 

In a first measurement, the spin states 4=(|1, — 1)±|1,+1)) 
were prepared and a small guiding field of 5.5 mG was ap- 
plied along the quantization axis z such that a slow oscillation 
(~ 8 kHz) could be observed. The remaining field compo- 
nents Bx and B, were compensated below 1 mG. After the 
prepared spin-states evolved a programmable period of time, 
the population of the ^(|1,— 1) — |1,+1)) state was mea- 
sured, see Fig. [3] (a). For such a field configuration, where 
B, dominates all other fields, the atomic state will stay within 
the subspace {|1, — 1),|1,+1)} during the Larmor precession, 
as |1,±1) are eigenstates of the interaction Hamiltonian. We 
observe the expected precession of an effective spin- 1/2 sys- 
tem with a 1/e dephasing time of about 120 /is. In order to 
extract the parameters responsible for the dephasing we have 
numerically fitted the dephasing model from Eq. ^ to the 
data points in Fig. [3] (a). This model includes fluctuations 
of the effective magnetic field consisting of residual fluctua- 
tions of external magnetic fields along the jc-axis (uniformly 
distributed) and the dominating distribution p{ABa,T) of the 
optically induced effective field ( fTSI l. For the fit we assumed a 
trap depth of Uq — 650 /iK and an average atomic temperature 
of 150/iK. The Larmor frequency deduced from this mea- 
surement corresponds to an average effective magnetic field 
component of B, = (5.5 ±0.5) mG composed of the optically 
induced field and the externally applied magnetic field. The 
observed dephasing is compatible with a standard deviation of 
the field distribution of 2.25 mG along the quantization axis. 
From this value we deduce a fraction of 0.6% of circularly 
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Figure 3: (Color online) Temporal evolution of different atomic 
spin states, a) Evolution of the superposition states 4=(|1, — 1) ± 
|1,+1)) in an effective magnetic field of 5.5 mG along the quanti- 
zation axis, b) Evolution of states |1,±1) in a field compensated 
to i? < 2mG. Measured were the populations of the spin states 
-T=(|l, — 1) — |l,+l))ina) and 1 1,-1-1) inb), respectively. The solid 
lines represent numerical fits of the measured data to the dephasing 
model in Eq. Q, which mainly incorporates fluctuations of the ef- 
fective magnetic field due to the vector light shift. 



polarized trapping light. This non-negligible fraction is due 
to the birefringence of the UHV glass cell where the experi- 
ment is performed. The birefringence is induced by mechani- 
cal stress and is not uniform over the walls of the cell, limiting 
the degree of control of light polarization at the position of the 
atomic trap. 



In a second measurement the evolution of the spin states 
|1,±1) was investigated. Here the absolute value of the effec- 
tive magnetic field was compensated to B < 2 mG. According 
to Fig. Ob), the stability of these states largely exceeds those 
of superpositions. This can be easily understood since the 
states |1,±1) are eigenstates of the effective magnetic field 
pointing along the quantization axis z, and therefore are not 
affected by the fluctuations along this axis. The slower de- 
phasing of the states |1,±1) shows that fluctuations along the 
X and y axes are smaller than along z. 
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Figure 4: (Color online) Partial tomographic reconstruction of the quantum state evolution. Shown are the density matrices (real part) for the 
prepared states 4=(|1, — 1)±|1, + 1)) in a) and b), and |l,Tl) in c) and d). 



E. Quantum state tomography 

The measurement of the temporal evolution of the atomic 
state provides a good way to determine its coherence proper- 
ties. However, state analysis in one basis does not give com- 
plete information about the qubit state under study. As the 
ground state 5^S\|2^F = 1 of ^^Rb has total angular momen- 
tum of one, the respective temporal evolution involves three 
Zeeman sublevels: nif — ±\ and nif = 0. Thus the analysis of 
dephasing processes becomes even more complex compared 
to a simple qubit state. 

The best way is a complete quantum state tomography (e.g. 
II28II and references therein), allowing to extract the informa- 
tion on how the state becomes mixed and also to distinguish 
coherent and incoherent processes. Unfortunately, a com- 
plete tomography of the spin-1 space in general requires 5 
Stern-Gerlach-like measurements (each providing the popu- 
lations of the 3 spin-1 eigenstates along a certain direction) 
II29I1 which are not accessible in our experiment. In particu- 
lar, the coherences between the 1 1 , ± 1 ) Zeeman states and the 
1 1 , 0) state can not be measured with our current technique, 
as the applied stimulated Raman adiabatic passage pulses an- 
alyze only the effective spin- 1/2 subspace {|1, — 1),|1,+1)} 
in a complete way 11191 ulll . However, as the detection effi- 
ciency is close to unity, we can infer the population poo of the 
1 1 , 0) state as the population missing in the {|1, — 1),|1,+1)} 
subspace. 

To reconstruct the density matrix p of the spin-1 ground 

ii2,F = 1, we use the worst-case assumption that 




there is no coherence between the |1,0) state and the |1,±1) 
states. We therefore set the corresponding off-diagonal den- 
sity matrix elements to zero. The corresponding full 3x3 
spin-1 density matrix is then given by 



(17) 



Here p^ is the 2x2 density matrix of the spin- 1/2 sub- 
space {|1, — 1),|1,+1)} and Poo = 1 —trace{p^^. As we typi- 
cally measure populations of the states |1,±1), -J=(|l, — 1)± 

1 1 , + 1 ) ) and 4=(|1, — 1)±;|1,+1)) which are eigenstates of 

Pauli spin operators ^z^^x^ and dy, the reduced density matrix 
Pv is fully accessible by our experimental techniques. By com- 
bining these complementary Stern-Gerlach measurements we 
are able to reconstruct the spin- 1/2 density matrix p, given by 



1 



I + {Ox) ■ Ox + (<7y) • 6y + (dz) • al 



(18) 



Now, in order to find a quantitative measure for the coherent 
fraction of the general spin-1 density matrix p , we decompose 
p as 



P - '■|Z)(Z| + (!-'■) 3 1, 



(19) 



state 5^S 



where |;t)(;t| is the density matrix of the closest pure state 
(which can be in general unknown), and ^1 represents a com- 
pletely mixed spin- 1 state. The corresponding pMnfy /?arame- 
ter r is the overlap with the closest pure state \x) and therefore 
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Figure 5: (Color online) Purity parameter resulting from the partial 
state tomography giving a lower bound for the state purity [6]. 



represents an ideal measure for the coherence of the investi- 
gated state. It can be calculated from the trace of p^ as 



'■=\l^i3-trace{p^)-l). 



(20) 



For r — Q the state under investigation is completely mixed, 
while for r = 1 it is a pure state. 

Based on the above procedure, tomographic measurements 
for the temporal evolution of spin-1 density matrices p for 
the initial states -!=(|1, — 1)±|1,+1)) and 1 1 , ± 1 ) were per- 
formed (Fig. m. The external magnetic field was set such that 
as little as possible Larmor precession could be observed up 
to 200 /is, and the circular fraction of the dipole light polariza- 
tion was < 1%. In the time evolution of the density matrices p 
of the initial states 4=(|1, — 1)±|1,+1)) one can observe sev- 
eral important features. The first one is the decay of the off- 
diagonal elements (coherences) as a general sign of dephas- 
ing. Second, a residual Larmor precession can be identified as 
the off-diagonal density matrix elements become imaginary 
(not shown in Fig. lU and undergo a change of sign. Third, 
the population of the 1 1 , 0) state continuously increases reach- 
ing ~ 15% after 200 /is, that is the qubit subspace is gradually 
depopulated. In contrast, for the initial states |1,±1) the ma- 
jor process during the evolution is only a slowly increasing 
population of the |1,0) state. 

In order to estimate the coherent fraction of the recon- 



structed spin-1 density matrices in Fig. |4]the corresponding 
purity parameter r was evaluated according to Eq. ( |20l i, as- 
sessing a lower bound of the atomic spin-coherence. For the 
evolution of the superposition states -U(|l, — 1)±|1,+1)), 

see Fig. |3a), we determine a 1/e dephasing time of 150/is. 
For the states |1,±1) we estimate the longitudinal dephas- 
ing time in absence of a guiding field by extrapolation to 
Ti > 500 /is (see Fig. |3b)). This value gives the time scale on 
which the populations of the effective spin- 1/2 states |1, — 1) 
and |1,+1) approach an equal mixture of all three spin-1 
basis-states |1,-1), |1,0) and |1,+1). 

In essence, the significantly longer longitudinal dephasing 
time shows that the fluctuations of the effective magnetic field 
are mainly along the quantization axis z- These fluctuations 
arise predominantly from the slow thermal motion of the atom 
in the trap where a residual circular polarization admixture 
leads to a position-dependent vector light-shift. The result- 
ing dephasing leads to a decay of the off-diagonal elements 
of the effective spin- 1/2 density matrix p^ with a 1/e time- 
constant of Tj = 75 /IS. The drift into the |1,0) state due to 
magnetic fields orthogonal to the quantization axis is signifi- 
cantly slower. 



IV. CONCLUSION 

In this paper we have studied the coherence properties of 
a qubit encoded in Zeeman substates of the hyperfine ground 
level of a single trapped **^Rb atom. While the "fundamen- 
tal" decoherence by Raman scattering of the photons from the 
dipole trap beam is negligible, the atomic state can dephase 
due to technical limitations. The main mechanisms leading 
to dephasing were identified as the fluctuations of stray mag- 
netic fields and the effective magnetic field induced by the 
circularly polarized component of the trapping light. By ana- 
lyzing the motion of the atom in the trap we have deduced the 
relation between the atomic temperature and the fluctuation 
of the effective magnetic field due to the circular admixture in 
the polarization of the trapping light. The dephasing of atomic 
memory states was then minimized by active stabilization of 
the external magnetic field together with an accurate setting 
of the polarization of the dipole trap light. 

By performing a partial state tomography of the 
5'^Si/2,F = 1 hyperfine ground level we have analyzed the 
dephasing of different states. The superposition states like 
-4=(|1,-1) ± |1,+1)) show a dephasing time of 75. .150/is 
which is mainly limited by field fluctuations along the quan- 
tization axis. The spin states |1,±1) are not sensitive to 
these fluctuations and thus show significantly longer dephas- 
ing times. Here we want to stress again that these dephasing 
times were measured at a magnetic field close to zero. An ex- 
ternally applied guiding field would induce a controlled pre- 
cession of the state while suppressing the influence of fluctu- 
ations orthogonal to its axis. However, the lifting of the de- 
generacy of the atomic states coming along with such guiding 
field may reduce the fidelity of the atom-photon entanglement 
scheme. Additionally it also would require a synchronization 
of the experiment (in particular of the time period between 



9 



preparation and measurement of atomic states) to the preces- 
sion period. 

In order to further extend the coherence time, two ways for 
improvement can be envisaged. On the one hand, better sta- 
bihty of the magnetic field can be reached by enlarging the 
geometry of stabilization coils, thereby reducing field gradi- 
ents. These measures may be combined with passive magnetic 
shielding for better suppression of external magnetic fields. 
On the other hand, a large contribution to the dephasing of the 
atomic ground state 5^5i m, F = 1 results from thermal motion 
of the trapped atom in the state-dependent potential induced 
by the residual fraction of circularly polarized dipole trap light 
(< 1%). Here longer coherence times could be reached with 
higher accuracy of the polarization alignment of the dipole- 



trap light and a reduction of birefringende of the glass cell, 
lowering of the trap depth, and/or better cooling of the trapped 
atom. Improvements of such type will extend the coherence 
times and thus the usability of neutral atom quantum memo- 
ries for future quantum repeater networks. 
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